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ABSTRACT 



The theorem of least work was used to investigate 
the stress pattern and bowing, i.e. uniform curvature, 
of an initially straight thin tube subjected to a steady 
temperature distribution independent of the longitudinal 
coordinate of the tube. Simple equations were developed 
for stress, strain, and bowing. They show that if the 
thermal strain is expressed as a Fourier series^ the tan- 
gential stress is small and is a function of the first 
harmonic only, the axial stress is a function only of 
the second and higher harmonics, and the bowing is a 
function only of the first harmonic . 
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1. Nomenclature. 



A = cross section area of tube wall 



a,b = coefficients used to specify axial strain 
(See Equation 4) 



£(«)- 



local thermal strain, presumed a function of e 
only, not of r or z. 



£*= coefficients in Fourier expansion of 



E = Young's modulus of elasticity 



n = an integer denoting the nth term of a 
Fourier series 



r = radial coordinate, measured positive outward 
from mid-thickness of tube 



ij, = tube radius to mid-thickness of wall 
t = tube wall thickness 

T(6) = local temperature 

T B = base temperature from which thermal strain is 
measured 

U = strain energy per unit length of tube 

u(0) = local radial displacement of mid-thickness 
curve of tube 

u„,u n = coefficients in Fourier expansion of u(0) 

v(0) = local tangential displacement of mid-thickness 
curve of tube 

v n ,v M = coefficients of Fourier expansion of v(0) 
x,y,z = tube coordinates (see Figure 1.) 

** = coefficient of linear thermal expansion 

ft = ^ + fVr^ 

^ elastic strain components in tangential and 

longitudinal directions, respectively 
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& - I + 4 r‘ 



= total strain (elastic plus thermal) in the 
tangential and longitudinal directions 



respectively 

t + 



n 

9 = 



angular coordinate measured clockwise from 
the positive z-axis, (see Figure 1). 



angular position of the plane of bowing 
K = curvature of bowed tube 
v = Poisson's ratio 



£ = /- J~ 



P n = / 



IZ, 






direct stress in the radial, tangential, and 
longitudinal directions, respectively 



shear stress in the r#~plane 

i+vCi-ftt ) 



Trd 

i. 

$ -fli + vO-V] 

x 

t 



*(!- Vll rj 

tt n t E 



_ /- i/‘ 






ui = 



_ d tv 

cfe 



TJ 



where oy is any representative variable 
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2. Introduction 

The problem of stress distribution in a tube with 
temperature variation independent of the axial coordinate 
has been solved under the following two sets of circum- 
stances: (1) an exact solution [2]* in the sense of the 
theory of elasticity for steady-state temperature config- 
uration, and (2) a very elementary soiution[5j in the 
framework of elementary strength of materials . Solution 
(1) is of limited usefulness, not only because of its 
limitation to the steady-state configuration but also 
because of its mathematical complexity., Solution ( 2 ) was 
developed to investigate technically important cases of 
the bowing of pipes partially filled with a cryogenic 
fluid, such as liquid oxygen. This solution assumes 
that only longitudinal strains and stresses in the pipe 
are an important cause of bowing . This in effect is 
taking Poisson’s ratio equal to zero. Presuming that 
Poisson’s ratio is zero, while actually it is in the 
neighborhood of 0.3 for most piping materials, led to a 
simple and useful solution for bowing radius and axial 
stress but said nothing about tangential stress. The 
degree of approximation involved was quite unknown; there- 
fore, it was the purpose of the present investigation to 
obtain a more accurate solution than that of Flieder, 
Loria, and Snith [53 without involving the mathematical com- 

^Numbers in brackets refer to works listed in the 
bibliography, page 45* 
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Then strain energy is differentiated partially with 
respect to each such coefficient, each such partial deri- 
vative is equated to zero, and finally the resulting set 
of equations is solved for the corresponding values of 
the coefficients. 

This procedure, which was used for the present investi- 
gation, gives the displacement which minimizes the strain 
energy (subject to the implied constraints of the assumed 
displacement) and thus corresponds to equilibrium. 

That the theorem of least work can be applied to the 
problem at hand, where temperature changes are involved, 
is not obvious. Argyris[l] and Boley and Weiner [2] have 
discussed, in general terms, the energy theorems and thermal 
stresses where temperature changes and non-linear elasti- 
city apply. Even so it seems appropriate to establish, 
by using the following physical considerations, that the 
theory of least work can be applied to the present problem. 

Consider the structure cut into a number of elemental 
volumes such that the temperature change is uniform in 
each element. A change in temperature causes thermal 
strains but not stresses; however, the elements no longer 
fit together continuously. In order to fit the elements 
together into a continuous structure having an assumed dis- 
placement pattern from the original shape, stresses must 
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be applied which cause elastic strains in accordance with 
the theory of elasticity. Since the stresses are applied 
holding temperature constant, the strain energy can be 
calculated using the equation previously cited, f<rd~r 
The values of the elastic strains are determined by the 
difference between the total strain, determined from 
geometric considerations of the assumed displacement pat- 
tern, and the thermal strains resulting from known tempera- 
ture changes. The theory of least work now applies and 
states that the strain energy is a minimum when the assumed 
displacement pattern is that corresponding to equilibrium. 
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3. Analysis of Stress and Strain. 

In order to proceed with the analysis, certain condi- 
tions must be met and assumptions made,, The tube must be 
constructed of an isotropic elastic material and the 
temperature distribution must be a known function of the 
angular coordinate, 6 , and independent of time and the 
longitudinal coordinate, z. This is the resulting tempera- 
ture distribution when a tube is partially filled with a 
cold liquid. Also no external forces or moments may be 
applied to the tube. The necessary assumptions are as 
follows: 

1. The tube is ,, thin n and thus and T r g are small 
and can safely be neglected in the analysis. 

2. The temperature is constant across the tube wall 
thickness. 

3. Young’s modulus of elasticity and Poisson’s ratio 
are constant. This is a reasonable assumption as shown by 
values for 302 stainless steel in Figures 1 and 2. 

Stainless steel is typical of most metals. 

4. Plane normal cross sections remain plane. 



By neglecting <T r and equilibrium cannot be satisfied 
and any results are necessarily approximate. A similar 
assumption has been made by others under similar circumstances 
and the solutions obtained were acceptable when compared 
with the exact solutions according to the theory of elasti- 
city. See DenHartog[4J pages 221 to 223 and page 240. The 
latter portion of Appendix III shows the relative magnitudes 
which may be expected of <r r and r r g . 
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5» In accordance with Saint Venant’s principle the 
stress pattern exists uniformly along the tube except at 
short distances from the ends. 




-400 -300 - ZOO -too O /OO 200 300 400 



TEMPER ATURE/f 



Figure 1 „ Modulus of Elasticity of AISI 302 
Stainless Steel 




Figure 2. Poisson 4 s Ratio for AISI 302 
Stainless Steel 



From T. F. Durham, R„ M. McClintock, and R„ P„ 
Reed, ’’Cryogenic Materials Data Handbook", 
National Bureau of Standards, PB Report 171809, 
1961 . 
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A. Symmetric Case. 

Temporarily the temperature distribution will be 
limited to one symmetric about the y-axis. The coordinate 
system used is shown in Figure 3. 



When there is a change in temperature of the tube, a 
point P located on the mid-thickness curve will, in gen- 
eral, be displaced to a new location Pi The radial dis- 
placement, u, and the tangential displacement, v, are 
positive in the directions indicated in Figure 4° 




x 



Figure 3 » Typical Section of Tube 




Figure 4» Displacement of Mid-thickness Curve 
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The tangential strain at any point can be expressed as the 
strain of the mid-thickness curve plus the strain due to 
bending of the tube wall, which in turn can be expressed 
in terms of the radial and tangential displacements as 

9 r ~ r&~ - (1) 



follows 

f o 



For derivation see Appendix I, 

Since the temperature distribution is symmetric about 
0 = 0 , the radial displacement will be an even function of 
$ and the tangential displacement will be an odd function 
of 0 . Expressed as Fourier series they are 

op 

u- Y co s (n&) 



As o 
oo 



( 2 a) 



( 2 b) 



and ( n9 '> . 

H = / 

The necessary derivatives are 

u = -Z nZu n cos(ne) 

V - Y n ^ cos(n&) . 

* - j 

Substituting Equations 2 into Equation 1 we obtain 

j: ! + + n IK,] cos (* &) . ( 3 ) 

The assumption that plane cross sections remain plane 
requires that the axial strain be 



= a + h(r 0 + r) cost'd) ) 



( 4 ) 



where a and b are constants to be determined, 
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Because there are no external forces or moments applied 



to the tube dA — 0 ( 5 ) 

£ y < r * dA=0 ( 6 ) 

f A *<*i dA=0 (7) 

where y = C r* + r) cos 6 and d A * + r) dr d Q . (g) 



Equation 7 is automatically satisfied as a consequence 
of symmetry. 

With the assumption that radial stress is zero, the 
stress condition in the tube is that of plane stress. 

The equations of plane stress are 

<r e = [et+ve'-O+^SM] (9) 

yf?d e * + v€ *~°*‘') £< ' e) ] , <io) 



where £(&)is the thermal strain. Expressed as a Fourier 

OQ> 

series g(e) = 8 a + 2 c °s (n9) , (11) 



where £, 



Z TT 



iir ,. T 

oz(t) <J~r 




and 



£ n - foe Cr)dT cos (nB)dB • 

o T-g 



The thermal strain is an even function as a consequence of 

the temperature being symmetric about the y-axis. In 
r T(«) 

general E(d)~ J 0 Using the above expression 

t b 
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for axial stress and integrating Equations 5 and 6 permits 
determining a and b in terms of displacement : 

« = <7^K>£ 0 - vJji* (12) 

Q 

b-fi[a +»)£,- (13) 

where £ = / ~ j^*. ,/S* / + and £ = / + ^ . 

The integration is performed in Appendix II . 

The elastic strain energy for the case of plane stress 
as discussed on page 4 is 

u = if{L € *-£te)]ok + [ee-£(6)]<fe} c/ >4 (H) 

per unit length of tube. 

Substituting Equations 9 and 10 for stress into Equation 
4, multiplying and collecting terms, we get 

= Jig+et-nve^-z-O+vlECSiese,] +zO+»)£(#f]jA. U5> 

Integrating Equation 15, we get the following expression: 

[z? ^ + n2 J 

+ i *r‘[2f4(io+/?b(“i+' / t )] iE 0 +fi& b£ t ] 

- 4f? &)£*] +zO+v)\z£.+'Ls5 ]' (16) 

See Appendix II for details of the integration. 
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Applying the theorem of least work as discussed pre- 
viously, we are now in a position to evaluate the coeffi- 
cients u ft and v n by taking the partial derivative of 'X 
with respect to each coefficient and equating the result 
to zero. 

For the coefficient u„, 







~4(!+i>)£jc-%£ ) ^ O „ 



Therefore, ^ (// o r 0 S 0 



where = — — £ . 

r * /->*<$ 



(17) 

(Id) 



For coefficients u, and v, , 

X- r 0 £b Z + -£*.&,+/,)* + &J?pb(u,+ v;) -Z6+1/)/6 k b£ t 

- zG+v) (u,+ i/j£, 
c 



where only those terir^s that contain u f or v 7 have been 
written., Note that b contains the sum ( u 7 +v # ) . It is 

not possible to evaluate u f and v, separately because 



equals c)X 0 However, it is sufficient to obtain their 
* u * dV t 

sum because that is the only form in which u, and v, appear 
as can be seen from Equation 19 . Taking and equating 
to zero: 



ML- .-i*Lb+ Vi&JL + )[i- -XaMi [ !-af\=o . , 20) 
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Solving for (u,+v, ) we get 

i 

(“,+ K) - ft, r 0 S, where 

For all other coefficients, u„ and v n , n>l : 

U y rt -r^v)^>,r 0 ^= o % 



( 21 ) 



and ttl/ -(H-^roS^ 



0 



yields n *4, 0+V)r o E n ~ /°n &n . 



( 22 ) 

(23) 

(24) 



Substituting Equation 24 into 22 yields (? n 'fn)U n =0 . 

But C < Pn~fn ) ^ 0 , therefore, 

; 14, « (l±vlp2Jk!L t (25a, b) 

The radial and tangential displacements are now 
known as functions of the thermal strain, and thus tempera- 
ture. Substituting Equations IS and 21 into Equations 12 
and 13 respectively, yields: 



cl- 



•— where £ — 1 + vO f (p Q ) (26) 

O 

and b = j*, Ej where “ -^-[f + v0~$)] m (27) 






Thus the axial strain is 



e 2 “ $ 0 £. +0+i )$,£,c°s(e). 



( 28 ) 
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Substituting Equations 18, 21, and 25 into 3 gives the 
tangential strain, 

ao 

o + 8 , c *s + 0+ k)X 6 n cos(ne) , ( 29 ) 

>7-2 

Remembering that £&> « £ 0 + E.cosCe ) +J* E n cos(ne) 

*~2j 

0O_ 

then 2 < 5 * cos (nO) = 8 (6) ~5 0 S, cos 8 and, 

ft-Z. 

e t ’‘Qtr')-v]£.+M-i)-v]e,cas6 +0*y)£(e). ( 30 ) 

Substituting these expressions! or strain into the stress 
formulas, Equations 9 and 10, determines the stress: 

< 5 ^= -E£( a)+ ^jt(£-0 ( 3 1 ) 

and <% = + »($>,- l)]S 0 -i-Hfl-t) +»($r0+»$, jz]S, cos &J , ( 32 ) 

This concludes the stress analysis for the symmetric case,, 

B. Anti -Symmetric Case* 

If the temperature distribution is anti “symmetric 
with respect to 0=0, the thermal strain is anti-symmetric „ 
Expressed as a Fourier series 

OO 

6 Ce) 8 n stn Cne) ( 33 ) 
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where 




Note that there is no constant term in the Fourier series. 

It is seen that this is similar to the symmetric case 
rotated 90 degrees and we can apply some of the results of 
the previous case. 

The first harmonic is the only one contributing to axial 
strain, so the axis of rotation is the y-axis and the 
longitudinal strain can be written as 



= a, + b, (r 0 -+rjs/n 6 > (34) 

where a, is automatically zero because of anti-symmetry. 

The radial displacement will be an odd function and the 
tangential displacement an even function of the angular 
coordinate. Expressed as Fourier series they are 

OO 

a.(9) =*£_ &„s/n (no) (35) 

/ 

oo 

and v(e)='E. v n cos(ne) . (36) 

Note that v contains no constant term. A constant term 
would imply a rotation of the tube. 

Substituting Equations 35 and 26 into 3» we obtain 



€^ = -£j[[[i+£( n - 0 ]u„ - Vu'jsmCne) . ( 37 ) 
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Letting ~ V n — ^ n yields 

\ +nX n }sm(nd) . OS) 

The expressions f or ^ and are of the same form as 

in the symmetric case except that sin(r?0) is substituted 
for cos{n9) « The only integrals that do not vanish when 
integrated over the interval 0 to 2tt are Js//j Z (nQ)d & <, 

These integrals have the same values as fco,s x (nO)d@ for 
the interval. Thus by making the substitution and 

writing si n(n&) for cos(«0), the solution of the symmetric 
case is applicable and the results are as follows: 



Q+n )$,£, sin 6 (39) 

e e = Sin 9 + (l+V)E (0) (40) 

c r z = -£E(d) + +v(f t -/) + £ ]s t sin d (4D 

<r 6 = jfp: fa, ~ 0 + +»$, -£]i,sin o , ( 42 ) 



The stress analysis is now complete for the anti-symmetric 
case . 
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C- General Case- 

For a general temperature distribution the thermal 
strain can be expressed as the sum of the symmetric and 



anti-symmetric cases, an even and an odd function- 

°° O© 

Thus £(&) = € 0 +^£„cos(n6) + ]E € n sm (n&) ) (43) 

a complete Fourier series- The general solution is as 
follows: € i =<& 0 £ o +0+j:)$j£ j cosO-i-S l s/n@] (44) 

€ e = fc&'O-vjBa +[(<£, -/) ->][£, cos 8 o] +(/+») e(e) (45) 

(%-~E E(6) + ~jjr{l(§ 9 ~f) + + 0-^)]6 o 

+ [£$,-0 4 vfy, ~0*0-^ *$,f][£/ cas & +£ l s/o&]}’ (46) 



°» = 7^^A-()+V^-/)j£ 0 +^ + l<? 1 -j)+l^{)[« | CM»4E 1 S/l.«j}(47) 

The stress and strain resulting from a general temperature 
distribution is now known- 

D- Simplified Equations for Stress and Strain- 

It is instructive to consider the magnitude of , 

(p o , § t , and lp t in order to simplify the equations of 
stress and strain- The stress and strain equations (Equa- 
tions 45 > 46 and 47) contain the terms {§ 0 -/ ), (^-/ ), 

(<^-/ ), and (^-/ ) which vanish if the ? s and ^ ? s are 
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unity o Table I shows that for tubes of practical wall 
thickness this is very nearly so. 

TABLE I 



The Effect of Wall Thickness on Stress and Strain 





0.1 


0.05 


0.01 


j". 


0.99992 


0.99996 


1.00000 




1 . 00027 


1.00012 


1.00005 


io 


1.00027 


1.00014 


1.00001 


$, 


0.99825 


0.99913 


0.99972 



for v = 0.296 



Substituting unity for ^ , $> f , and 



Equations 44 , 45, 46, and 47 become: 

€ z = 6 o +0+-£ )(S,cos9 +€,S/ri & ) (43) 

e d = 6+»)£(e) -p(E* +£,cos e +£, s //'? <9 J ( 49 ) 

OJ = -e\e(Q) + (50) 

0" =^T ^ fey cos £ + £, 9) . (51) 



These equations should be sufficiently accurate for most 
design applications . 
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4. Bowing of Tube. 



A tube subjected to the assumed temperature distribu- 
tion will bow into the arc of a circle, or if restrained 
develop stresses equal to those required to straighten the 
bowed tube. This is of practical concern in pipe design 
and has been discussed by Flieder, Smith, and Wetmore[6j. 
From the analysis of the bending of an initially straight 
beam [7] it is known that the curvature, K , is^- ; where 
c is the distance from the neutral axis and € is the strain 
resulting from bending. The axial strain on the mid- thick- 
ness curve of the tube is &+£,&/» 6) ? obtained 

by taking r=0 in Equation 44 . The term denotes a 
general elongation of the tube and does not contribute to 
bowing. The S, c °£ & term results in bowing in the yz-= 
plane. The strain, € , at ® *0 is $>£>, at a distance r 0 

from the neutral axis. Thus , Ku t ~ is the curvature 

in the yz-plane. The £> t 6.,s/n9 term results in bowing in 
the xz-plane. The strain at is and as before the 

distance to the neutral axis is ri . Therefore, Kxt s 



is the curvature in the xz-plane. The curvature in general, 



is the vector sum of and Ayz-, and the plane of bowing 




(52) 



is given by d H - , 



(53) 
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measured clockwise from the y-axis. Again it is possible 
to set =1 without significant error so that 







r 0 



( 54 ) 
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5. Results, Design Equations. 

The analysis leads to the following simplified equa- 
tions which can be used to design thin tubing systems 
subjected to steady temperature distributions that are in- 
dependent of tube length. 

1. Axial strain 



€ a = + )(E,cos& +E,s/n&) 



(4*) 



2 . Tangential strain 



€ ff s 6+v)£(6) -v(e o +s,cos0 + E,sme) 



(49) 



3. Axial stress 




4. Tangential stress 



ft cos0 + 8,s/na) 



(51) 



5 . Bowing curvature 




(54)&(53) 








rvc TC&J 



Where 



z7r S T ($J 



£ m = oc(r)dT cos(nO)(Jg 



° T s 



_2X T(Q) 



E«~Z7rf f <*( r )dT sm(n&) dd 



° t b 
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6. Discussion of Results. 

When thermal strain is expanded in a Fourier series, 
the tangential stress and bowing as given by Equations 51 
and 54 are functions of the first harmonic only, and the 
axial stress as given by Equation 50 is a function of the 
second and higher harmonics. To interpret the physical 
meaning of this, consider the case where £(o ) = £, 

When the thermal strain is projected onto a yz-plane,the 
first harmonic transforms to a linear function and the 
second harmonic transforms to an even function. It is 
known [8] that linear thermal strain does not cause stress 
in a free beam, but causes bowing, and that an even thermal 
strain function causes stress but not bowing. Thus, the 
results are those expected from superimposing elementary 
solutions for simple beams. 

Equations 44 thru 47 include a first approximation for 
the effect of wall thickness; however, the basic assumptions 
of this analysis, particularly that the temperature does 
not vary through the wall, suggest that the results be 
confined to thin walled tubes. Accordingly, Equations 4& 
thru 51 and 54 are suggested as being the appropriate ones 
to use for actual design and analysis. These results 
afford a slight extension of those obtained by Flieder, 
Loria, and Smith in the following respects: (1) they give 
a variation of axial stress through the wall thickness and 
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(2) they give a value for tangential stress. 

The present results and those of Flieder, Loria, and 
Smith are comparable for the case where E, V , and oc are 
constant. When these conditions are imposed, Equation 50 
is identical at the wall mid-thickness to the result devel- 
oped by Flieder, Loria, and Smith. Note that 0^ is 
independent of Poisson’s ratio. Equation 50 gives a minor 
variation of through the wall; they did not consider 
any radial variation of . Flieder, Loria, and Smith do 
not give a result for tangential stress. Equation 51 shows 
that tangential stress is zero at the tube wall mid-thickness 
and that it is small for tubes with slight bowing. Equation 
51 also shows that tangential stress vanishes if Poisson’s 
ratio is zero. 

Comparison of the results for bowing show that there 
is a very slight effect due to finite wall thickness and 
non-zero Poisson’s ratio as given by Equation 52. However, 
the simplified equation, Equation 54» which is here recom- 
mended for purposes of design and analysis, is essentially 
identical to that given by Flieder, Loria, and Smith. 

Essentially, we have shown that influences neglected 
in the analysis of Flieder, Loria, and Smith have only the 
very slightest effect on the final result. However, the 
expression of the thermal strain as a Fourier series has 
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simplified the final form of the equations and aids in the 
visualization of the effect that a given temperature dis- 
tribution may be expected to have on a tube. 
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simplified the final form of the equations and aids in the 
visualization of the effect that a given temperature dis- 
tribution may be expected to have on a tube. 
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APPENDIX I 



Derivation of Expression for Tangential Strain 

The temperature change from point to point around the 
circumference of the tube results in thermal strains. There 
is assumed to be no temperature gradient across the wall 
thickness of the tube, thus the thermal strain is indepen- 
dent of radius. Superimposed upon this strain is another 
strain which varies with radius and is caused by the change 
in curvature of the tube wall. 

To derive an expression for that portion of the strain 
independent of radius, consider two points A and P, close 
together, on the mid-thickness curve. When the temperature 
changes, the points will move to new locations A, and P ( as 
shown in Figure 5. 




x 



Figure 5. Tangential Strain at Mid-thickness. 
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Let /• and 7 be unit vectors in the i and j directions 



respectively. The bar denotes a vector. 

Then OA, = vl + (r Q +u)j j (la-1) 

oq = (r 0 +tz+do)\sm(dQ)l + cos(d&)j] > (la-2) 

QP, = (v+d v)[cos (c/0) 7 + s/n(d6)j] ? ( la _ 3 ) 

and A,f) ~ OQ QFf ~~ OA f > (la-4) 

Because dd is a small angle, sin (dd) =d@ and cos (dd) =1. 
Then, Atf -fa+u.)de+ dud&-t-dvjj +ldu-vd9- dvd&Jj . (la-5) 

Neglecting infinitesimals of higher order, 

A, P, ~ [(Z+u)d&+dv]T +[du + vdO) J (la-6) 

and \A,?d=L( r *+ u ) X d^+zfo+<i)vdd Z +CL*'dO l -Zu.vdO i +v\/0J . (la-7) 



The strain, , on the mid-thickness curve is defined 
as the change in length per unit initial length of the curve 
as the tube wall deflects. For small angles the arc length 
and the cord length can be taken as equal. Therefore, the 



strain is 



e e,/ 



_ M,?f - rs c/<9 

r„d& 



(la-8) 



Substituting the value of l A ,Vl yields 






_ fC&ufd&i +zln«i) irdd^+ ULdB^- ZuvA V^J - f 0 d& 

r 9 de 



(la-9) 



Again neglecting infinitesimals of higher order, the strain 
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becomes r<> [/ + ~ (u.+ v) ] e — f 0 




(la-10) 



O 




(la-11) 



Therefore , 



- VkL v 
'&>! K 



(la-12) 



on the mid-thickness curve. 

The second component of strain is caused by the change 
in curvature of the tube wall* Because the tube is "thin" 
the theory for bending of straight beams can be used* The 
strain due to bending is 

F &A , (la-13 

where r = distance from mid-thickness and 

= change in curvature of tube wall mid-thickness 



To find the change in curvature consider a point P 
which after a change in temperature is at a new location 
P, as shown in Figure 6* 



curve o 
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y 




Figure 6. Displacement of a Point on Mid-thickness 
Curve . 

The point is located by the vector R= (C> +CL )er+ V&e (la-14) 

where & r and are unit vectors in the radial and tan- 
gential directions respectively. Taking derivatives; 

# = ue r + (r a +u) e* + v G-e ~ v e r 

= (u-v)e r +(rz + u+ v)e a (la-15) 

and S - (a-v)e r +(ul- V)e& + (u+ V)e B - fa+u. 4 1 >) e r 

= C(X-Zv-0L-r o )e r +(V+Zu~v)ee . (la-16) 
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The curvature after deflection isfeJ 



i/_ L £L -Zv-U‘-KXr* + u.+ v) M 1rn 

• ( 71 

Noting that r e greatly exceeds u, v and their derivatives, 
we get the approximation 

‘ ^ jH. 

= r o z -( 0 L+&)r o . (la-18) 

r 3 

f o 



The change in curvature is the final curvature minus the 
initial curvature. The tube is initially round so the 
initial curvature is l/r 0 . Thus, 



,K= r /-( u +“) - 1 = - 



r* 



6 






(la-19) 



The strain due to bending is then 




(g±u ) r 

C z 



(la-20) 



The tangential strain is the sum of the two components 
developed above : 



• - \L _ (U-+ il C 
■9 r 0 r 0 *• 



(l) 
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Solving Equation 2a-4 for a we have 



a. - 0^)8 o - L 1 ( 2a-5) 



and <Z- 0+ y)£ 0 ~~ 


r 0 


7 . 

£ = / . , ^ ■*. 
where / /< r « . 


(12) 


Equation 6. 








£<%ydA=0 






(6) 


y = (r 0 +r)cos 9 


and 


dA=(fo+ r ) drd@ 


(S) 



Substituting Equations 8 and 2a-l into 6 and integrating 
over the same interval as before: 

— drj[cLcos0+ b(r 0 + ^cos^O +VZ.T [L i +z( n ' , )fao+ * 
l-» i 

~C)-hi02iSjzo^C^O)cosO]clO = 0, (2a-6) 

> 1*0 ^ 



f** r tt * * 

Note that J cosc*e) cos(n9)d0 . 



(2a-7) 



Henceforth in the analysis, terras which vanish in this way 
will not be written. Continuing the integration, 

S - ? 2 

I [b (r 0 +rf +y(u,+ V)(z±jj 2 - (/+»)&, (r 0 +r) 1 ’]dr = 0 (2a-8) 

' 4 * 
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Let <f- / 



= / + 






tz c* > 



a rs °- h< * ~ 

then b = ~A~ [0+v)£, - vCu.. + K) ] 

<* £ r a 



Equation 15 . 

zo-pv = f ~f\+el -&+»)£&) ^ 

-2-(]+v)8(d)€ 9 + zO+v)e(efl JA 



( 13 ) 



(15) 



The integral is the sum of the following terms: 

I / e* JA = fr.-tr)j r J[a}+ bVr„+r/^ s ^j 

^ -to 

x 

= fWc) fe 71 " &+ if b x (r a + r)^]d r 

'i 

= £~TT GL lr Q {+ 7 TbV» 3 ^ f/ * ] 

* zira^ni -hir^-tlz £ 

A *7T 

^/} sr fe+r)drf^^{[l 4 "£/» x '!)] «£ +Znd + '6’>-d]u n ^ + *T t&JcAStfisWd 



IZT 



'4 



o ° o 



r % 

= J (&*)?* fa -£)£+ Z +2"C'+ ?&-')]** 



-K 



ZZi f u? +Z u«-t2nAU-*it,+ *K)\ 

C I 
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I 



where £ = I+ /Ti; z ^~^ a>act /°r, = ,+ • 

k r i7r 

€ & JA = ^ C(r^r)JrJ[Ci-^)a.a 0 + Cr o -i-r)Cu^,) bcos x 0lct0 



/a. 



zm: 

- v r J 

o ■' 

= SJjJT J (r 0 +>-)lz(i+ %)au.', + fr.+r)(u l +v,)b]dr 

r * y ~k 

- ZyzLIsJ L £ + (u,+ v>) bj 8 

^ Jit 

JV ' f-Z{/-hv)£(d)e i clA= -26+y)f (r b +r)Jrfl<3.£s+ b(K+ r )8 l cos z 0]d 6 

\ 



= -2 Oi-y)yr fiz a.^+ 6 Cr a +*j£, ](V*r) ^ r 
'A 

= -Z0+v)7Tr o t [zaSo + b Toft £,] 



*/ z ** 

f-Z6+»)S(e) e 9 dA = ~ 2 ^ > J ( r >+ r )drf^$+ffr-')K + * &}z£os x (h6)d# 

A ° ° 

oo 

<o -it/ 1 



= -zty&Si &U.S. +t(f.Uj-» K)£j 



£ *r 

JZL J zQ+»)Ste) Z dA = 2 0+y)f.fc+r)d^y el cos * 6 *®) d-e 

A ~K~o 



= zO-n>)-*Q t [2 e * +Z e*„ ] 
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Summing all terms: 



ZO-W - 

7Tg i£ 



=Za z + + ^,[2 £u? (^u* +2 np n V n +n'K \ ? 

+Zr[2GUo f+Cu? v)bfi]-2(i + »)L&-£ a + b/3r 0 6,] 



-z^[z$Uv£ 0 +^(fiu.»+nVjSn J 



+zO+v)[2£o +1^£«] . 



( 16 ) 
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APPENDIX III 



Example 

• The following example will illustrate the method of 
solution and order of magnitude of results. A 302 stain- 
less steel pipe 10 inches in diameter, with a wall thick- 
ness of 0.5 inches is subjected to the temperature 
distribution shown in Figure 3. Assuming the pipe is free 
from restraint, what is the stress pattern and how much 
will the pipe bow? 




e, RADIANS 



Figure 3. Temperature of Pipe. 
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The first step is to determine the thermal strain for 
each temperature . The base temperature in this case is 
taken as 6£°F because the data, see Figure 9, is given in 
that form. A harmonic analysis is then performed to deter- 
mine £o > S, j and S, . The calculations and results are 
given in Table II and paragraphs below. 




TEMPERATURE , °F 

Figure 9. Thermal Expansion of AISI 302 Stainless Steel 

From T. F. Durham, R. M. McClintock, and R. P. Reed, 
"Cryogenic Materials Data Handbook", National 
Bureau of Standards, PB Report 171$09, 1961. 




TABLE II 



Harmonic Analysis 



8 


T(0) 


SCO) 


CO 

0 

0 

§ 


£6^3 sin 


radians 


°F 


X/0* 


x/o^ 


x/o r 


0.0 


50 


-20 


-20.0 


0.0 


0.1 


43 


-21 


-20.0 


- 6.5 


0.2 


43 


-26 


-21.1 


- 15.3 


0.3 


35 


*35 


-20.6 


-23.3 


0.4 


23 


-45 


-14.2 


-43.3 


0.5 


-10 


-62 


0.0 


-62.0 


0.6 


-20 


-30 


24.7 


-76.2 


0.7 


-70 


-125 


73.5 


- 101. 4 


0.3 


-250 


-260 


211.0 


*153.0 


0.9 


-295 


-235 


272.0 


-33.2 


1.0 


-300 


-233 


233.0 


0.0 


1.1 


-295 


>235 


272.0 


33.2 


1.2 


-250 


-260 


211.0 


153.0 


1.3 


-115 


-162 


95.3 


131.2 


1.4 


-100 


-150 


46.3 


142.5 


1.5 


-35 


-140 


0.0 


140.0 


1.6 


-75 


-130 


- 40.2 


123.3 


1.7 


-50 


-110 


- 64.7 


39.2 


1.3 


15 


-50 


- 40.5 


29.4 


1.9 


40 


-30 


-23.5 


9.3 




Total 


-2566 


1224.0 


323.0 



!.9V 

E 0 = ^ II €(e) = - IZ 8- 3 Xfo' 

/•d'tf 

E, « TslE(e)cose = izZA x /a 

o 

£, = 7^ <? (&Js//7 e=*3Z.3x /o' 



The bowing of the pipe is 



— Z.^3x/0 /In 



n> 



at 



O IP — /a.*-' S, - i&n 1 32A =. /4. 8 . 



IZ2-4 
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The tangential stress is given by 



$e ^ e sm 9 ^ > 

therefore, on the mid-thickness curve, r=0, <7g- 0. 

Letting E=30xl0 6 psi and >'=0.296, 
vA. _ o.z9i(3o*) _ q 7 ^ 6 

/-^ ,-o.o&iz ~ y ' /6 e Sj - • 

At the outside diameter r=t/2=0.25 inches. 

C T q - 9.7 cos# + 32-3 s/v 0 ) t 

The results are given in Table III. 

The axial stress is given by 

-E{s(&) ^ ,+ iS^vy Jfe °ose -t- i, sw9) ]} 

= -[-iz&f+fi + r tfrP))(l2Z'fcos d +3Z-ifsme) ]} # 

On the mid-thickness curve r=0 and [ I ^ = 1» 

At the outside diameter r=0„25 and [/+ 3 - 1 - 055 - 

Axial stress was calculated for each position and the 
results are given in Table III and also shown in Figure 10. 
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TABLE III 



■ 



Tangential and Axial Stress 



6 


,psi 


, psi 


0i , psi 


radians 


at r=0.25 


at r=0 


at r=0.25 


0.0 


596 


4,200 


5,700 


0.2 


575 


4,800 


8,700 


0.4 


333 


“4,200 


-3,000 


0.6 


-23 


”16,500 


-16,800 


0.8 


“390 


14,400 


13,500 


1.0 


-596 


11,400 


9,300 


1.2 


”575 


4,200 


2,400 


1.4 


”333 


-13,800 


-15,000 


1.6 


23 


2,700 


3,000 


1.8 


m 


600 


2,700 




Figure 10. Stress Pattern of Pipe. 
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The result of the first assumption, that <T r and T r $ 
are negligible, is that equilibrium cannot be satisfied. 
The values of these stresses are now estimated using the 
conditions of the example. A section is cut out as shown 
in the free body diagram, Figure 11. These are the posi- 
tions of maximum and zero tangential stress. 



y 




Figure 11. Shear Stress 

The stress distribution, , causes a net counter-clock- 

wise moment of 25.7 in-lbs. To restore equilibrium, 
stresses 0^ and 7? re are required. Summing moments about 
the point P, 

7 ro*Mih ■ /0 - 2 P s ‘ ■ 

The tangential stress, <r Q ' , is also equal to 10.3 psi. 

An indication of the magnitude of the radial stress 
required to maintain equilibrium can be obtained by inves- 
tigating the condition at 0=195°, the cold side of the pipe. 
There the axial stress is about 10,000 psi and will be 
taken as constant. Figure 12 shows the action of axial and 
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radial stress. 




<%= 10,000 pS( 



Figure 12. Radial Stress Component Caused by 
Axial Stress. 



The radius of curvature, R, was obtained from the 
bowing. Summing forces in the radial direction, 

‘R(T r d6- {(T^dd ~0 . 

Therefore, <V ' ^ § ' ^ 0 °°° ^ = 1.27 psi ’ the radial 

stress caused by axial stress. 

The maximum radial stress component resulting from 
tangential stress can be estimated by considering the stress 
condition at S- dp as shown in Figure 13. 
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Figure 13 . Radial Stress Caused by Tangential 
Stress. 

Taking the sum of the forces in the radial direction 
we get, 

z<r r r 0 J&- a ^C\<r d + %' )6d 
^ r ~ °i / § j + /0 ) ~ 776 P sc • 

The total radial stress is the sum of the component 
caused by axial stress and the component caused by tangen- 
tial stress ( G~r ~ •= 9 ps < ) . 

It is seen from the above calculations that op. and 

are small when compared with other stresses, and there- 
fore, neglecting them is realistic. 
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